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Abstract 

The uniqueness of the orthogonal Z 7 -circle patterns as studied by Bobenko 
and Agafonov is shown, given the combinatorics and some boundary conditions. 
Furthermore we study (infinite) rhombic embeddings in the plane which are qua- 
sicrystallic, that is they have only finitely many different edge directions. Bicoloring 
the vertices of the rhombi and adding circles with centers at vertices of one of the 
colors and radius equal to the edge length leads to isoradial quasicrystallic circle 
patterns. We prove for a large class of such circle patterns which cover the whole 
plane that they are uniquely determined up to affine transformations by the com- 
binatorics and the intersection angles. Combining these two results, we obtain the 
rigidity of large classes of quasicrystallic Z 7 -circle patterns. 

1 Introduction 

Circles, especially circle packings and circle patterns, have successfully been used over 
the past years to define and study discrete analogs of classical smooth objects. In 
particular, this approach leads to discrete holomorphic mappings, for example discrete 
analogs of the power functions z 1 , and to discrete holomorphic function theory. See for 
example [231 El HI IH] for some of the contributions to the theory of circle patterns and 
[2"T] for results on circle packings. 

In this article we focus on circle patterns which are characterized by a given combi- 
natorics specifying which circles should intersect and by the corresponding intersection 
angles. Thus we associate to a circle pattern a pattern of kites corresponding to in- 
tersecting circles, see Figures [T] (right) and [6] A particularly suitable source for the 
required knowlegde on circle patterns and their relations to consistency, integrability, 
and discrete holomorphic functions is the textbook in discrete differential geometry. 
Our main results can roughly be summarized as follows: The combinatorics (together 
with suitable intersection angles and boundary conditions, if necessary) determines the 
geometry of the circle pattern. This rigidity result can be interpreted as a discrete 
version of Liouville's Theorem in complex analysis. 

Our first result concerns the uniqueness of circle patterns which are discrete analogs 
of the power functions z 1 for 7 6 (0,2) as defined in [71 01 [2] • Here, the square grid 
combinatorics, the orthogonal intersection angles, and the two boundary lines of a sector 
uniquely determine the geometry of the pattern, see Figure [2] for an illustration. The 
rigidity of orthogonal Z 7 -circle patterns was only known for rational 7, see [2]. 

Furthermore, we consider the case of a circle pattern which covers the whole complex 
plane and for which the radii of all circles are equal and the interiors of different kites 
are disjoint. Then the corresponding kites form a rhombic embedding. Also, assume 
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that there are only finitely many different edge directions of the kites. Such rhombic 
embeddings are called quasicrystallic [5] • Additionally, orient an edge e and consider 
a line perpendicular to the edge. Move this line parallelly in positive and in negative 
direction along e. We suppose that in both cases this moving line intersects infinitely 
often edges parallel to e. We assume that this property is true at least for two edges with 
linearly independent directions. Such rhombic embeddings of the plane are for example 
given by Penrose tilings, see for example Figure [6] We show that any other embedded 
circle pattern with the same combinatorics and intersection angles of the kites is the 
image of this embedding by an affine transformation. 

Rigidity of some classes of infinite circle patterns of the plane have already been 
studied. Schramm |23j considers square grid combinatorics and orthogonal intersection 
angles. As an essential step of our proof we generalize his result to square grid circle 
patterns with regular intersection angles ip £ (0,7r) and (tt — ip). He [19] studies disk 
triangulation graphs and exterior intersection angles in [7r/2,7r] which does not cover 
the class of isoradial quasicrystallic circle patterns defined above. 

As observed in [5j isoradial quasicrystallic rhombic embeddings can be used to define 
corresponding quasicrystallic Z 7 -circle patterns, see also Examples of quasicrystal- 
lic Z 7 -circle patterns are shown in Figures [7J |9j and [Tl] They have been created using 
software developed by Veronika Schreiber for her diploma thesis |24j . Our rigidity result 
for orthogonal Z 7 -circle patterns is also generalized for large classes of quasicrystallic 
Z 7 -circle patterns. 

There is some more literature concerning rigidity for infinite planar circle packings, 
that is configurations with non-overlapping touching circles, see [5TJ [TSJ [32 |5D] . Our 
rigidity proofs adapt some of the ideas which have been used for packings. In particular, 
we apply discrete potential theory. 

This paper is organized as follows. First we introduce terminology and present useful 
facts about circle patterns. We especially focus on regular circle patterns with square 
grid combinatorics. Then we recall in Section[4]the definition and some properties of the 
orthogonal Z 7 -circle patterns for 7 S (0, 2) as studied in [H[TJ[2] and prove their rigidity. 
In Section [5] we introduce quasicrystallic circle patterns and prove uniqueness for a class 
of these patterns. Finally, we recall the definition and some facts about quasicrystallic 
Z 7 -circle patterns for 7 e (0, 2) as studied in [3J [§1 HJ and prove rigidity for certain 
classes of these patterns. A more detailed version of the results can be found in [12") . 

2 Circle patterns 

In this section we focus on a definition and some useful properties of circle patterns. We 
describe circle patterns using combinatorial data and intersection angles. 

The combinatorics are specified by a b-quad-graph Si, that is a strongly regular cell 
decomposition of a domain in C possibly with boundary such that all 2-cells (faces) are 
embedded and counterclockwise oriented. Furthermore all faces of Si are quadrilaterals, 
that is there are exactly four edges incident to each face, and the 1-skeleton of S is a 
bipartite graph. We always assume that the vertices of Q> are colored white and black. 
To these two sets of vertices we associate two planar graphs G and G* as follows. The 
vertices V(G) are all white vertices of V(%>). The edges E(G) correspond to faces of @, 
that is two vertices of G are connected by an edge if and only if they are incident to the 
same face. The dual graph G* is constructed analogously by taking as vertices V{G*) 
all black vertices of 2$. Si is called simply connected if it is the cell decomposition of a 
simply connected domain of C and if every closed chain of faces is null homotopic in Si. 

For the intersection angles, we use a labelling a : F(S) — > (0, n) of the faces of S. 
By abuse of notation, a can also be understood as a function defined on E(G) or on 
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Figure 1: Left: An example of a b-quad-graph @ (black edges and bicolored vertices) 
and its associated graph G (dashed edges and white vertices). Right: The exterior 
intersection angle a of two intersecting circles and the associated kite built from centers 
and intersection points. (3 = ir — a is the interior intersection angle. 

E(G*). The labelling a is called admissible if it satisfies the following condition at all 
interior black vertices v € Vi nt {G*): 



Definition 2.1. Let @ be a b-quad-graph with associated graph G and let a : E(G) — > 
(0, 7r) be an admissible labelling. An (immersed planar) circle pattern for @ (or G) 
and a are an indexed collection ^ = {C z : z 6 V(G)} of circles in C and an indexed 
collection X = {K e : e £ E(G)} = {K f : f e F(&)} of closed kites, which all carry the 
same orientation, such that the following conditions hold. 

(1) If z 1 ,z 2 G V(G) are incident vertices in G, the corresponding circles C Zl ,C Z2 in- 
tersect with exterior intersection angle a([zi, Z2]). Furthermore, the kite K[ Zl ,z 2 ] ^ s 
bounded by the centers of the circles C Zll C Z2 , the two intersection points, and the 
corresponding edges, as in Figure [l] (right) . The intersection points are associated 
to black vertices of V($$) or to vertices of V(G*). 

(2) If two faces are incident in the corresponding kites share a common edge. 

(3) Let fi, . . . , f n E F(3)) be the faces incident to an interior vertex v € Vi n t(3$). Then 
the kites K* x , . . . , Kf n have mutually disjoint interiors. The union Kf 1 U • • • U Kt n 
is homeomorphic to a closed disk and contains the point p(v) corresponding to v in 
its interior. 

The circle pattern is called embedded if all kites of have mutually disjoint interiors. 
The circle pattern is called isoradial if all circles of ^ have the same radius. 

Note that we associate a circle pattern ^ to an immersion of the kite pattern J(f 
corresponding to @ where the edges incident to the same white vertex are of equal 
length. The kites may also be non-convex and can be constructed from a (suitable) 
given set of circles and from the combinatorics of G. 

As all circle patterns considered in this article will be planar and immersed, the 
notion "circle pattern" will include these properties in the following. 

For our study of a circle pattern c € we will use the radius function r<# = r which 
assigns to every vertex z € V(G) the radius r<g{z) — r(z) of the corresponding circle 
C z . The index ^ will be dropped whenever there is no confusion likely. The following 
proposition specifies a condition for a radius function to originate from a planar circle 
pattern, see [TU] for a proof. 





(1) 



/ incident to v 
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Proposition 2.2. Let G be associated to a b- quad- graph S> and let a be an admissible 
labelling. 

Suppose that c € is a planar circle pattern for 2> and a with radius function r = r<^ . 
Then for every interior vertex Zq G Vint (G) we have 

( E /a([*,*o])a°gr(z)-logr(z ))J "T = 0, (2) 
\[z,z ]eE(G) I 

where 

1 1 — e x ~~ i0 

and the branch of the logarithm is chosen such that < fg{x) < tt. 

Conversely, suppose that & is simply connected and that r : V(G) — > (0, oo) satis- 
fies ([2]) for every z e Vi nt (G). Then there is a planar circle pattern for G and a with 
radius function r. This pattern is unique up to isometrics of C 

For the special case of orthogonal circle patterns with the combinatorics of the square 
grid, there are also other characterizations, see for example [23] . 

Note that 2/ Q ,([ ZZo ])(logr(z) — logr(zo)) is the angle at zq of the kite with edge 
lengths r(z) and r(z ) and angle a([z,Zo]), as in Figure [I] (right). Equation ^ is the 
closing condition for the closed chain of kites which correspond to the edges incident to 
Zq. This corresponds to condition (3) of Definition |2.1| 

For further use we mention some properties of fe, see [261 Lemma 2.2]. 

Lemma 2.3. (i) The derivative of fg is f' g (x) = 2(cos h'"- cos g) > 0. 
(ii) The function fg satisfies the functional equation fg{x) + fg(—x) = n — 9. 

If there exists an isoradial circle pattern, we can obtain another circle pattern from 
a given radius function. 

Lemma 2.4. Let G be a graph constructed from a b- quad- graph 2> and let a be an 
admissible labelling. Suppose that there exists an isoradial circle pattern for G and a. 
Let r be the radius function of a planar circle pattern for S> and a. Then there is a 
circle pattern ^ for G and a with radius function r<g = 1/r. 



Proof. By Lemma 2.3 (ii) the function 1/r satisfies condition ([2| of Proposition 2.2 for 



all interior vertices z$ € Vi n t(G). In particular, 

1 , 1 

Mo]6£(G) 



Y f a ([z,z ])(tog-rx - lo g-r^) 
L — ' u " r{z) t(Zq) 



Y - a([z,z })) - Y /«(M)( lo K^j) = 2 ' r -' r = ' r - 

[z,z ]£E(G) [z,z ]eE(G) ^ 0> 

Here we have also used that ^ ]g£;(G)( 7r — a ([ z i z o])) = 27r since there is an isoradial 
circle pattern for G and a and the assumption that r is the radius function of a circle 
pattern for G and a. □ 

Let G be a graph constructed from a b-quad-graph @ and let a be an admissible 
labelling. Suppose that ^1 and ^2 are planar circle patterns for @ and a with radius 
functions n = r<g\ and ri = tv 2 respectively. Define a comparison function w : V(S>) —¥ 
Cby 

iw(y) = r 2 (y)/r 1 (y) for y G V(G), 

\w(x)=e lS ^ eS 1 ioixeV{G*). U 
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Here S(x) £ R or w(x) = e lS ^ is defined to be the rotation angle or the rotation 
respectively of the edge-star at x € V(G*) when changing from the circle pattern c €\ to 
^2- Note that w(y) is the scaling factor of the circle corresponding to y £ V(G). Then 
w satisfies the following Hirota Equation for all faces / £ F(&). 

iv(x )w(y )a - w(xi)w(y )ai - ti;(a;i)w(yi)ao + w(x )w(y 1 )ai = (4) 

Here xo,xi £ V(G*) and yo,yi £ V(G) are the black and white vertices incident to 
/ and ao — xo — yo and a\ = x\ — ya are the directed edges. Thus equation Q is 
the closing condition for the kite corresponding to the face /. Furthermore, the Hirota 
Equation is 3D-consistent; see Sections 10 and 11 of [9] or [11] for more details. This 
property will be used in Section [5j 



3 SG-circle patterns 

In this paper we are particularly interested in the special case of regular circle patterns 
with square grid combinatorics. First, we fix some notation. Let SGD be the regular 
square grid cell decomposition of the complex plane, that is the vertices are V(SGD) = 
Z + ill and the edges are given by pairs of vertices [z, z'\ with z, z' £ V(SGD) and 
\z-z'\ = 1. The 2-cells are squares {z + a + ib : a,b £ [0,1]} for z £ V(SGD). As SGD 
is a b-quad-graph, the vertices 

V(SG) = {n + im£Z + iZ:n + m = (mod 2)} 

are assumed to be colored white. As above, SG and its dual SG* are defined as the 
associated graphs to SGD. Furthermore, SG(n,v) with nfi and v £ V(SG) denotes 
the subgraph of all vertices with combinatorial distance at most n from v in SG. 

Let ip £ (0,7r) be a fixed angle. Define the following regular labelling on E(SG). 
Let [z\,Z2] be an edge connecting the vertices z\, Z2 £ V(SG). Without loss of generality, 
we assume that Re^) < Re(z 2 ). Then 

a ([zi z 2 \) = if ImOi) < Im(z 2 ), ^ 

' l 7r — V' if Im(z!) > Im(z 2 ). 

If G is a subgraph of SG, a circle pattern for G and is called SG-circle pattern. 
The choice tp = n/2 leads to orthogonal S'G-circle patterns as considered by Schramm 
in [23]. 

Theorem 3.1 (Rigidity of S'G-circle patterns). Suppose that c & is an embedded planar 
circle pattern for SG and . Then ^£ is the image of a regular isoradial circle pattern 
for SG and under a similarity. 

The proof is a suitable adaption of the corresponding proof for orthogonal S'G-circle 
patterns given by Schramm using suitable Mobius invariants, see [23] Theorem 7.1] 
or [T2J. This adaption needs the following generalization of the Ring Lemma of [2T] . 
which is also useful in the following. 

Lemma 3.2. Let r be the radius function of an embedded circle pattern for SG(3, 0) 
and a^. There is a constant C = C(ip) > 0, independent ofr, such that for k = 0, 1,2,3 
there holds 

r(0) > 
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Proof. Assume the contrary. Then there is a sequence of embedded circle patterns 
for SG(3,0) and such that r„(0) = 1 and r n (i k (l + i)) — > as n — > oo for some 
k 6 {0,1,2,3}. Without loss of generality we assume that k — 0. We also may as- 
sume that the circle Co corresponding to the vertex G V(SG) and the intersection 
point corresponding to 1 € V(SG*) are fixed for the whole sequence. Then there is a 
subsequence such that all the circles converge to circles or lines, that is converge in the 
Riemann sphere C = S 2 . Now equation ^ implies that there exist some kites which 
intersect in the limit in their interiors. But this is a contradition to the embeddedness 
of the sequence. □ 

If the number of surrounding generations is big enough, there is the following useful 
estimation on the quotient of radii of incident vertices. 

Lemma 3.3. There is an absolute constant C > such that the following holds. 

Let G be a subgraph of SG and let ^ be an embedded circle pattern for G and 
with radius function r. Let v € V(G) be a vertex and suppose that SG(n,v) C G, that 
is G contains n generations of SG around v, for some n > 3. Then for all vertices w 
incident to v there holds 

C r(w) , C 

l--<-Vf<l + -. 6 
n r[v) n 

The proof is an adaption of the corresponding proof for hexagonal circle packings 



given by Aharonov in and uses Lemma 3.2 see [T^] for more details. The necessary 



results on discrete potential theory can be found in the appendix of 13J and in [12 



4 Uniqueness of orthogonal Z 7 -circle patterns 

An orthogonal circle pattern with the combinatorics of the square grid associated to the 
map z 1 was introduced by Bobenko in j? j . Further development of the theory is due to 
Agafonov and Bobenko [H [TJ [2] . 



4.1 Definition and useful properties 

In the following we briefly summarize the definition and some known facts about or- 
thogonal Z 7 -circle patterns, see [H [I] [5] for more details. 

Definition 4.1. Let fl C 2 2 . A map / : D —> C is called discrete conformal if all its 
elementary quadrilaterals are conformal squares, i.e. their cross-ratios are equal to —1: 

Q(,fn,m : fn-\-l,m: fn+l.ra+1 1 fn,7n+l ) ■ 

(/n,m /n+l,m ) (/n+l,m+l /n,m+l) 

(,/n+l.m ./n+l.m+l ) {fn,m+l fn,m*) 

Here and below we abbreviate f n , m = f{n, m). 

A discrete conformal map f rhm is called embedded if the interiors of different elemen- 
tary quadrilaterals (/„, m , /„+i, m , /n+i,m+i, fn,m+i) are disjoint. 

Note that the definition of a discrete conformal map is Mobius invariant and is 
motivated by the following characterization for smooth mappings: 

A smooth map / : C C D — > C is called conformal (holomorphic or antiholomorphic) if 
and only if for all z = x + iy G D there holds 

lim q(f(x,y), f{x + e,y), f(x + e,y + e), f(x,y + e)) = -1. 



G 



Figure 2: Illustration of the discrete conformal map Z 3 / 2 and the orthogonal Z 3 / 2 - 
circle pattern (right). 



In order to construct an embedded discrete analog of z 1 the following approach is 
used. Equation ([T]) can be supplemented with the nonautonomous constraint 

r ~ (/n+l,m fn,m)yfn,m fn — l,m) 

Tjn,m — ^ n 7~7 7 \ 

\Jn+l,m jn—l,m) 

. n \fn.m+\ fn.mjifn.m /n.m— l) /0 \ 

+ 2m — — . (8) 

V7n,m+1 Jn,m—1) 

This constraint, as well as its compatibility with Q, is derived from some monodromy 
problem; see @]. We assume that < 7 < 2 and denote 

Z 2 + = {(n, m) EZ 2 :n,m> 0}. 

The asymptotics of the constraint Q for n, m — > 00 and the properties z 7 (lR + ) = 
R+ and z~*(iM+) = e~< m / 2 R + of the holomorphic mapping z 1 motivate the following 
definition of the discrete analog. 

Definition 4.2. For < 7 < 2, the discrete conformal map Z 1 : 1? + — > C is the solution 
of equations Q and pi with the initial conditions 

Z 7 (0,0)=0, ZT(1,0) = 1, Z 7 (0,1) =e 77ri/2 . 

From this definition, the properties Z^(n, 0) €E M+ and Z 7 (0,m) € e im / 2 M+ are ob- 



vious for all n,m G I. Furthermore, the discrete conformal map Z 1 from Definition 4.2 



determines an S'G-circle pattern. Indeed, by Proposition 1 of [4] all edges at the vertex 
fn.m with 7i + to = (mod 2) have the same length and all angles between neighboring 
edges at the vertex f n . m with n + m = 1 (mod 2) are equal to 7r/2. Thus, all elemen- 
tary quadrilaterals (/„ >m , f n +i,m, /n+i.m+i, /n,m+i) build orthogonal kites and for any 
(77, to) € Z\ with 77 + 77?, = (mod 2) the points f n +i,m, fn,m+i, f n -i,m, fn,m-\ lie on a 
circle with center f n , m - Therefore, we consider the sublattice {(n, to) 6 1? + : n + m = 
(mod 2)} and denote by V the quadrant 

V = {z = N + iM : N,M 6Z,M> |iV|} C Z 2 , 

where N = (n — to)/2, M = (77 + to)/2. Two vertices Zj, ^ 6 V are connected by an 
edge if and only if \zi — z^\ = 1. 

Theorem 4.3 ([H[IJ[2]). (i) If R(z) denotes the radius function corresponding to the 
discrete conformal map Z 1 for some < 7 < 2, then it holds that 

(■y-l)(R(z) 2 -R(z-i)R(z + l))>0 (9) 

for allzeW\ {±N + iN\N € IN}. 
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(ii) For < 7 < 2, the discrete conformal maps Z 1 given by Definition 4-2 are embed- 



ded. Consequently, the corresponding circle patterns are also embedded. 

In the following section and in Section [6] we continue to use the notation of this 
section. In particular the radius function is denoted by R and we have the normalization 
R(0) = 1. 

4.2 Uniqueness of the orthogonal Z 7 -circle patterns 

This section is devoted to the proof of following uniqueness result. 

Theorem 4.4 (Rigidity of orthogonal Z 7 -circle patterns). For 7 6 (0, 2) the infinite or- 
thogonal embedded circle pattern corresponding to Z 1 is the unique embedded orthogonal 
circle pattern (up to global scaling) with the following two properties. 

(i) The union of the corresponding kites of the Z 1 -circle pattern covers the infinite 
sector {z = pe i/3 £ C : p > 0, (3 £ [0, 771-/2]} with angle jir/2. 

(ii) The centers of the boundary circles lie on the boundary half lines R+ and e l77r / 2 lR,+ . 

Our proof uses results of discrete potential theory or of the theory of random walks 
which can be found in standard textbooks, for example by Doyle and Snell ,14] or by 
Woess [28] . We recall some basic terminology and notation and cite adapted versions of 
a few theorems which will be useful for our argumentation. 

By abuse of notation, we denote by 1? the points (a, b) £ R 2 = C with a, b £ Z as 
well as the graph with vertices at these points and edges e = [zi,z 2 ] if | ^1 — 1 = 1. 
The meaning will be clear from the context. 

Consider the network (Z 2 ,c) with conductances c(e) > and resistances l/c(e) on 
the undirected edges e £ Eifl?). Then a transition probability function p is given by 

, „ x ._ Sc([zi,z 2 ])/(j: e=[zuz]eE{z2) c(e)) if [z u z 2 ] £ E(Z 2 ) 
I otherwise 

The stochastic process on 1? given by this probability function p is a reversible random 
walk or a reversible Markov chain on 1? . The simple random walk on 1? is given by 
specifying c(e) = 1 for all edges which leads to p(z\, z 2 ) = 1/4 if [zi,z 2 ] £ E(fl?~). 

Denote by p e sc the probability that a random walk starting at any point will never 
return to this point. The network (Z 2 ,c) is called recurrent if p csc = (and transient 
otherwise). Note that p csc = 1/R e g, where R e g denotes the effective resistance from a 
point to infinity. 

Theorem 4.5. (i) The simple random walk on 1? is recurrent. 

(ii) Let (Z 2 , Ci) and (Z , c 2 ) be two networks with conductances Ci(e) > and c 2 (e) > 
on the edges. If c 2 (e) < c\(e) for all edges e £ E(fE 2 ), then the recurrence of 
(Z 2 ,ci) implies the recurrence of(1i 2 ,c 2 ). 

A proof can for example be found in 14, Chapters 5, 7, and 8] or [23 Sections l.A, 
l.B, and Corollary (2.14)]. 

A function / : Z 2 — > R is called superharmonic (subharmonic) with respect to the 
probability function p or with respect to the conductances c if for every vertex v £ Z 2 

we have J2 we ^ p( v > w )f( w ) < /O) (Etuez* P( v > w )f( w ) > f( v ))- 

The following proposition shows that the quotient of the radius functions of two 
orthogonal SG-circle patterns is subharmonic with respect to suitably chosen conduc- 
tances. As the statement is a special case of Proposition |6.5| below, we omit the proof. 
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Proposition 4.6. Consider two orthogonal circle patterns for SG(1,0). Denote the 
radii by pj and Tj respectively, where po and ro denote the radii of the inner circles. 
Then 

3 = 1 ' 3 j=l ^ 3 = 1 3 3 = 1 U 

where Cj = l/((pj/p ) + (po/pj))- 

Our proof of rigidity is based on the following property of superharmonic functions 
on recurrent networks. 

Theorem 4.7 ([28, Theorem (1.16)]). A network is recurrent if and only if all nonneg- 
ative superharmonic functions are constant. 

Proof of Theorem \4-4\ Let 7 <E (0,2) and denote by R : Y — > R + the radius function 
of the embedded Z 7 -circle pattern with i?(0) = 1. Let r : Y — > R + denote the radius 
function of an embedded orthogonal circle pattern with the same combinatorics and 
the same boundary conditions (orthogonal boundary circles to the half lines R+ and 
e* 77r / 2 R + ). Without loss of generality we assume same normalization r(0) = 1. This 
can always be achieved by a suitable scaling. In the following, we will show that the 
radius functions R and r take the same values on all of Y. This implies that both circle 
patterns coincide. 



As both circle patterns are embedded, Lemma 3.3 implies that for some constant 
A > and n > 3 



" r(*F) 



A r(z) n >) A , N 



holds for all radii r(zg ) f° r vertices Zq of the nth generation away from the origin and 

(n) 

their incident vertices Zj . Here, vertices z G Y belong to the nth generation if their 



combinatorial distance in Y to the origin is n. For estimation (111 we have also used 
that the reflection of the circle pattern in one of the boundary lines R + or e 477r / 2 R + also 
leads to an embedded orthogonal S'G-circle pattern. The same reasoning applies to the 
radii of the Z 7 -circle pattern, so 

A R{zf ] ) A , s 

1 - - < ] . < 1 + - 12 
" " R(z { n) ) n 

for n > 3 with the same constant A. Estimati ons pTj ) and ( 12 ), the boundary conditions 
and a suitable adaptation of the Ring Lemma |5.11| for circles of generation two and three 
from the origin imply that there is a constant K > such that 

K r(z Q ) K R(zo) 

for all incident vertices Zj and zq. 

We now consider two undirected networks (Z 2 ,C) and (Z 2 ,C) as follows. On the 
edges of Y we define two conductance functions C and C by 



C(e) = C(R( Zj ),R(z k )) :-- 
C(e) = C(r(z j ),r(z k )):= 



(R( Zj ) R(z k ) 
\R{z k ) + R( Zj ) 

r ( z 3) , r ( z k) 

r{z k ) r(zj) 
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where the edge e = [zj, Zk] connects the vertices Zj, Zk € V. Estimations ( 13 ) imply that 
both positive functions C > and C > arc uniformly bounded away from (and from 
infinity). These two conductance networks on V C Z 2 can be continued to all of 1? by 



reflection in the lines {\M\ = \N\}. From Theorem 4.5 we deduce that both networks 
(Z 2 ,C) and (Z 2 ,(7) are reccurent. 

Consider the following positive functions on Y 

f 1 (z) = r(z)/R(z)>0 and f 2 {z) = R{z)/r{z) = l/f x {z) > 0. 



By Proposition |4.6| these functions are subharmonic. Using the boundary conditions 
of the circle patterns, this remains true if fi and f 2 are continued to all of 1? using 
reflection. Consequently, M — f\ and M — f 2 are superharmonic for all constants Met. 
If /i or f 2 is bounded from above, we get a positive superharmonic function using the 



upper bound. Then Theorem 4.7 implies that both functions are constant. Thus r = R 
and consequently both circle patterns coincide. 

Denote by Mi(n) and M 2 (ri) the maximum of fi and f 2 , respectively, for the set 
of vertices of the nth generation about the origin. As fi and f 2 are subharmonic, 
they assume their maxima on the boundary. Therefore the functions Mi and M 2 are 



monotonically increasing. The estimations (11) and (12) imply that the quotients of 
any two radii of one circle pattern in the nth generation are bounded from above for 
n > 3, as two vertices in the nth generation can be connected by at most 4n edges using 
only vertices of the nth and n + 1st generation. So their quotient is bounded by e AA for 
both radius functions r and R. Note that with the normalization r(0) = 1 = R(0), the 
maxima Mi and M 2 are bounded from below by 1. Thus their product 

(nh _ r ( z mI) r ( z< M 2 ) , m 



Mi{n)M 2 {n) = h(z$)f 2 (z$) = ^ 
Ml 2 r(*<g) R{z (£) 



(n)\ - 



is bounded from above. Here z^ 1 and z$ denote the vertices of the nth generation 
where /i and f 2 assume their maxima, respectively. Therefore Mi and M 2 are also 
bounded. This finishes the proof of uniqueness. □ 



5 Uniqueness of isoradial quasicrystallic circle pat- 
terns 

The uniqueness result of Theorem |3.1| can be generalized for some classes of quasicrys- 



tallic circle patterns. On this basis we will then generalize Theorem 4.4 for some classes 
of quasicrystallic Z 7 -circle pattern. 



5.1 Quasicrystallic circle patterns and connection to Z 

Definition 5.1. A rhombic embedding in C of a b-quad-graph Si is an embedding with 
the property that each face of S is mapped to a rhombus. Given a rhombic embedding 
of ^, consider for each directed edge e € E(&) the vector of its embedding as a complex 
number with length one. Half of the number of different values of these directions is 
called the dimension d of the rhombic embedding. If d is finite, the rhombic embedding 
is called quasicrystallic. 

Adding circles with centers in the white vertices of the rhombic embedding and 
radius equal to the edge length reveals the close connection to embedded isoradial circle 
patterns. 

A circle pattern for a b-quad-graph Si is called a quasicrystallic circle pattern if there 
exists a quasicrystallic rhombic embedding of S and if the intersection angles are taken 
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from this rhombic embedding. The comparison function of the isoradial circle pattern 
for and the quasicrystallic circle pattern ^2 is also called comparison function for 

c &2- 



In the following we will often identify the b-quad-graph $ with a rhombic embedding 



of 




Figure 3: An example of a 
quad-surface £l%> C 1? . 



Remark 5.2. The notion "quasicrystallic" is not uniquely defined in literature. Here we 
adopt the definition given in [S]. Naturally, this property only makes sense for infinite 
graphs or sequences of graphs with growing number of vertices and edges. 

Any rhombic embedding of a b-quad-graph 'S can 
be seen as a sort of projection of a certain two-dimen- 
sional subcomplex (quad- surface) Q,g> of the multi- 
dimensional lattice 7L d (or of a multi-dimensional lat- 
tice C which is isomorphic to 1i d ). An illustrating 
example is given in Figure [3] 

The quad-surface f2® in Z d can be constructed 
in the following way. Denote the set of the different 
edge directions of ^ by A = {±ai, . . . , ±<Zd} C S 1 . 
We suppose that d > 1 and that any two non-opposite 
elements of A are linearly independent over R. Let 
ei , . . . , e^ denote the standard orthonormal basis of 
R d . Fix a white vertex xq G V(3>) and the origin 
of WL d . Add the edges of {±ei, . . . , ie^} at the ori- 
gin which correspond to the edges of {±di, . . . , ±ad} 
incident to Xq in @, together with their endpoints. Successively continue the construc- 
tion at the new endpoints. Also, add two-dimensional facets (quadrilateral faces) of 1i d 
corresponding to the quadrilateral faces of f^, spanned by incident edges. 

A quad-surface f2@ in Z d corresponding to a quasicrystallic rhombic embedding can 
be characterized using the following monotonicity property. For a proof see 9, Section 6]. 

Lemma 5.3 (Monotonicity critcrium). Any two points of can be connected by a 
path in with all directed edges lying in one d- dimensional octant, that is all directed 
edges of this path are elements of one of the 2 d subsets of {±ei, . . . , ±e<j} containing d 
linearly independent vectors. 

An important class of examples of rhombic embeddings of b-quadgraphs can be 
constructed using ideas of the grid projection method for quasiperiodic tilings of the 
plane; see for example [16j |TTl [25] . 

Example 5.4 (Quasicrystallic rhombic embedding obtained from a plane). Let E be 
a two-dimensional plane in R d . Let ei,.. .,e<j denote the standard orthonormal basis 
of R d and let t G E. We assume that E does not contain any of the segments Sj — 
{t + Aej : A € [0, 1]} for j = 1, . . . , d. HE contains two different segments Sj 1 and 
Sj 2 , the following construction only leads to the standard square grid pattern Z 2 . If 
E contains exactly one segment Sj, the construction may be adapted for the remaining 
dimensions (excluding ej). We further assume that the orthogonal projections onto E of 
the two-dimensional facets Ej 1 ^ 2 = {Aie^ + A 2 e j2 : Ai, A 2 G [0, 1]} for 1 < ji < J2 < d 
are non-degenerate parallelograms. Then we can choose positive numbers C\, . . . , such 
that the orthogonal projections PE(cjej) have length 1. 

Consider around each vertex p of the lattice C = c\% x • • • x c^Z the hypercuboid 
V = [— ci/2, ci/2] x • • • x [— Crf/2, Cd/2], that is the Voronoi cell p + V. These translations 
of V cover R d . If E intersects the interior of the Voronoi cell of a lattice point (i.e. 
(p + V)° l~l E j£ for p G £), then this point belongs to tt c (E). Undirected edges 
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correspond to intersections of E with the interior of a (d— l)-dimensional facet bounding 
two Voronoi cells. Thus we get a connected graph in C. An intersection of E with the 
interior of a translated {d — 2)-dimcnsional facet of V corresponds to a rectangular two- 
dimensional face of the lattice. By construction, the orthogonal projection of this graph 
onto E results in a planar connected graph whose faces are all of even degree (= number 
of incident edges or of incident vertices) . A face of degree bigger than 4 corresponds to 
an intersection of E with the translation of a (d — fc)-dimensional facet of V for some 
fc > 3. Consider the vertices and edges of such a face and the corresponding points and 
edges in the lattice C. These points lie on a combinatorial fc-dimensional hypercuboid 
contained in C By construction, it is easy to see that there are two points of the fc- 
dimensional hypercuboid which are each incident to fc of the given vertices. Choose a 
point with least distance from E and add it to the surface. Adding edges to neighboring 
vertices splits the face of degree 2fc into fc faces of degree 4. 

Thus we obtain an infinite monotone two-dimensional quad-surface fl c (E) which 
projects to an infinite rhombic embedding covering the whole plane E. Parts of such 
rhombic embeddings are shown in Figures |6j [HJ and |10| 

5.2 Quasicrystallic circle patterns and integrability 

Let @ be a quasicrystallic rhombic embedding of a b-quad-graph. The quad-surface 0@ 
in Z d is important by its connection with integrability. See also [IT] for a more detailed 
presentation and a deepened study of integrability and consistency. 

In particular, a function defined on the vertices of which satisfies some 3D- 
consistent equation on all faces of can be uniquely extended to the brick 

n(n§>) := {n = (m, . . . , n d ) e Z d : a fc (ft@) < n k < b k (Q^), k = l,...,d}. 

where a fc (f2@) = min ne y(n @) n k and 6fc(Slg>) — max ng v (!)9 ) n k . Note that II(fi@) is the 
hull of 0®. A proof may be found in |9l Section 6]. Let 9> be a quasicrystallic rhombic 
embedding and let ^2 be a quasicrystallic circle pattern with the same combinatorics and 
the same intersection angles. Denote the comparison function for ^2 by w as in ^ . Since 
the Hirota equation Q is 3D-consistent (see Sections 10 and 11 of [5]) w considered as a 
function on V(£Lg) can uniquely be extended to the brick 11(51®) such that equation Q 
holds on all two-dimensional facets. Additionally, w and its extension are real valued 
on white points of V(d@) and have values in S 1 for black points of V(tt@). This can 
easily be deduced from the Hirota Equation Q. 

The extension of w can be used to define a radius function for any rhombic embedding 
with the same boundary faces as @. 

Lemma 5.5. Let &) and S>' be two simply connected finite rhombic embeddings of b- 
quad-graphs with the same edge directions. Assume that 3> and 'S' agree on all boundary 
faces. Let c € be an (embedded) planar circle pattern for S> and the labelling given by the 
rhombic embedding. Then there is an (embedded) planar circle pattern for & and 
the corresponding labelling which agrees with for all boundary circles. 

Proof. Consider the monotone quad-surfaces and f2^. Without loss of generality, we 
can assume that and Sl'g, have the same boundary faces in 1i d . Thus both define the 
same brick II(f2@) = n(il^) =: II. Given the circle pattern < ^', define the comparison 
function w for c if on V(O^) by Extend w to the brick LI such that condition Q 
holds for all two-dimensional facets. Consider w on fi^ and build the corresponding 
pattern c &' , such that the points on the boundary agree with those of the given circle 
pattern . Equation Q guarantees that all faces of fl'^ are mapped to closed kites. Due 
to the combinatorics, the chain of kites is closed around each vertex. Since the boundary 
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kites of c €' are given by which is an immersed circle pattern, at every interior white 
point the angles of the kites sum up to 27r. Thus c &' is an immersed circle pattern. 

Furthermore, is embedded if ^ is, because c iS' is an immersed circle pattern, and 
c €' and ^ agree for all boundary kites. □ 

5.3 Local changes of rhombic embeddings 

Let Si be a rhombic embedding of a finite simply connected b-quad-graph and let 0@ 
be the corresponding quad-surface in Z d . Let z € Vi n t(Q@) be an interior vertex with 
exactly three incident two-dimensional facets of Consider the three-dimensional 

cube with these boundary facets. Replace the three given facets with the three other 
two-dimensional facets of this cube. This procedure is called a flip; see Figure [4] for an 
illustration. 



Figure 4: A flip of a three-dimensional cube. Only the faces bounded by solid edges 
are part of the quad-surface in 1i d . 

A vertex z € Z d can be reached with flips from if z is contained in a quad-surface 
obtained from Q@ by a suitable sequence of flips. The set of all vertices which can be 
reached with flips, including V(rt@), will be denoted by J-"(f2g>). 

Remark 5.6. The quad-surface ft c (E) can also be generalized using a finite sequence 
of flips. Such an infinite rhombic embedding will be called a plane based quasicrystallic 
rhombic embedding. 

As a generalization of simple flips we define flips for simply oder doubly infinite strips 
of the following form. See Figure [5] for an illustration. 




Figure 5: An example of a flip for an infinite strip. 

Let fi® C Z d be a simply connected monotone quad-surface. Let z € V w (fl%t) be 
a white vertex. Let , e j2 , e j3 be three different edges incident to z such that there 
are two-dimensional faces /i,/2 of incident to e 3l and e j2 , and to ej 2 and e^ 3 , 
respectively. Let a\ = a(fi) and a-i = a(/2) be the intersection angles associated to 
these faces. Let «3 be the intersection angles associated to the two-dimensional facet 
of Z d incident to Bj 1 and ej 3 . Then X)j=i a i ~ ^ n or X^=i( 7r ~ a i) + a 3 = 27r. In 
the first case, consider the half-axis ^+ = {z + Xej 2 : A > 0}, where e j2 is the vector 
corresponding to the edge e j2 and pointing away from z as in Figure [5] In the second 
case, consider the other half-axis i^L = {z + Ae j2 : A < 0}. In both cases we may also 
consider the whole axis S% = {z + Ae^ : A € R}. Assume that the translations of f\ 
and /2 along these (half-)axis, that is the faces fj + nej 2 + z for j = 1, 2 and n G IN, 
(—n) £ IN, or n £ Z respectively, are contained in f2@. We only consider the case of 
the positive half-axis ^+ further. For the argumentation is analogous and the case 
of the whole axis Si is a simple consequence. Replace each face /i + nej 2 + z by its 
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translate f\ + ne.j 2 + z + §j 3 for n £ JNo and similarly fa + nej 2 + z by fi + n ^j 2 + z + 
for n G INo, where and ej 3 are the vectors corresponding to the edges and ej 3 
respectively and pointing away from z. Adding the face incident to z, e^, and ej 3 , we 
obtain a different, but still monotone simply connected quad-surface. 

The definition of an infinite flip for a black vertex z £ V;,(f2g>) is very similar. 

Lemma 5.7. Let fig> C Z d 6e a simply connected monotone quad-surface and let f2'^ be 
the simply connected monotone quad-surface obtained from f2@ after performing a flip 
for a simply infinite strip. Let ^ be a circle pattern for & and the corresponding labelling 
and let c £' be the corresponding circle pattern after performing the corresponding infinite 
flip as for Q'@ . Then the resulting circle pattern "jf is embedded if the original one < to is. 

The proof is based on similar arguments as the proof of Lemma |5.5| and is therefore 
left to the reader. 



5.4 Uniqueness of isoradial quasicrystallic circle patterns 

Let S be the family of all infinite quasicrystallic rhombic embeddings of connected and 
simply connected b-quad-graphs *3i which cover the entire complex plane and such that 
the brick LI(Sl^) of the corresponding quad-surface contains a Z 2 -sublattice, that 
is there are at least two different indices Ji,j2 such that min^y^^) n j k = — 00 an d 
max ngv(f! 9 ) n j k = 00 f° r k = 1,2. Note that <§ contains in particular the plane based 
rhombic embeddings, like the Penrose tilings, for which n(Sl^) = X d . Now we use the 
uniqueness of S'G-circle patterns of Theorem |3.1| in order to establish the uniqueness of 
the circle patterns of S '. 

Theorem 5.8 (Rigidity of quasicrystallic isoradial circle patterns). Let & £ § be an 
infinite quasicrystallic rhombic embedding with associated graph G and corresponding 
labelling a. Let be an embedded circle pattern for G and a. Then ^ is the image of 
the isoradial circle pattern corresponding to & under a similarity of the complex plane. 

Proof. Let S> € <§ be an infinite quasicrystallic rhombic embeddings with associated 
graph G and corresponding labelling a. Let ^ be an embedded circle pattern for G 
and a. Consider the comparison function w for ^ defined by ([3| on and extend 
it to LT(f2@). Let z £ V(f2@). By assumption on there is a Z 2 -sublattice f2(z) 
with z £ V(il(z)) which is contained in 11(0®). Furthermore, we can perform flips 
for fig> and corresponding flips for the circle pattern ^ such that the resulting quad- 
surface f2' contains an arbitrary number of generations of the lattice S!(z) about z. As 



the corresponding circle pattern c € l is embedded by Lemma 5.7 and as the number 



of generations about z can be chosen arbitrarily large, we deduce from the Rigidity 



Theorem 3.1 that the radius function is constant on Q(z). More precisely, the extension 
of w is constant on white and black vertices of f2(z) respectively. 

This argumentation is valid for all vertices z £ V(fl^) and all Z 2 -sublattice f2(z) with 
z £ V(ft(z)) which are contained in LT(f2@). Therefore the extension of w is constant, 
on white and black vertices respectively, on all Z 2 -sublattices which are contained in 
LT(r2@). Due to our assumptions on the combinatorics of 3, the radius function which 
is w restricted to white vertices has to be constant on the whole brick fl(Sl^). This 
implies in particular that is an isoradial circle pattern and thus is the image of the 
isoradial circle pattern corresponding to 3 under a similarity of the complex plane. □ 

Definition 5.9. Let 3 be a rhombic embedding of a simply connected b-quad-graph. 
The combinatorial K -environment of a point z £ V(£t) is the subgraph corresponding 
to all vertices which have combinatorial distance at most K from z in 
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Corollary 5.10. Let S> G S be an infinite quasicrystallic rhombic embeddings with 
corresponding labelling a. Let vq G V w {2$) be a white vertex. Then there are a constant 
no = no{2>) G IN and a sequence s n [vo,&) decreasing to for n — > oo such that the 
following holds. 

For n G IN, n > no, let 3$2n(vo) be the rhombic embedding corresponding to the 2n- 
environment of Vq. Let G n (vo) be the associated graph. Let %?„ be an embedded circle 
pattern for G n (vo) and the labelling a taken from 3> with radius function r n . Then there 
holds 

r n [vo) 



1 



<s n (v ,@) (14) 



r„(vi) 

for all vertices v\ G V{G n (vo)) incident to vo- 

We omit the proof which is very similar to the proof of the Hexagonal Packing 



Lemma of [21j . The following lemma is similar to Lemma 3.2 and corresponds to the 
Ring Lemma of [2"T] . 

Lemma 5.11. Let %> be a quasicrystallic rhombic embeddings with associated graph G 
and labelling a. Denote by a m i n = min{a(e) : e G E(G)} the smallest intersection angle. 
Let no G IN be such that (no — 3)a m i n > tt. Let vo G V W (S>) be a white vertex. Assume 
that & contains a (2no) -environment about Vq. Then there is a constant C = C(S>) > 
such that the following holds. 

Let r be the radius function of an embedded circle pattern for & and a and let V\ be 
a vertex incident to vq in G. Then 

r(vo) 

Proof. Suppose that there is a vertex v\ incident to vo and a sequence of embedded circle 
patterns for and a with radius functions r n such that r„(u ) = 1 and r n (vi) — > as 
n — > oo. Without loss of generality we may assume that the circle Co corresponding 
to the vertex vq and the intersection point corresponding to one fixed black vertex Wo 
incident to vo and v\ in & are fixed for the whole sequence. Then there is a subsequence 
such that all the circles converge to circles or lines, that is converge in the Riemann 
sphere C = S 2 . 

Equation ^ implies that there are at least two circles corresponding to vertices 
incident to vq whose radii do not converge to 0. If the limit is finite, there are at least 
two circles whose radii do not converge to corresponding to vertices incident to this 
vertex of the first generation and so on. 

Consider the kites which contain in the limit the intersection point corresponding to 
Wo and apply the above argument at most no times. Then by our assumption on n 
we obtain two kites whose interiors intersect in the limit. This is a contradition to the 
embeddedness of the sequence. □ 

If @ is a plane based quasicrystallic rhombic embedding we also obtain an analog of 
the Rodin-Sullivan Conjecture, see [2"T1 ITBl IS]. 

Corollary 5.12. Let Si be a plane based quasicrystallic rhombic embedding. There are 
absolute constants C = C(£>) > and no — no(S>) G IN, depending only on S>, such that 
for all white vertices vo G V W (S>) and all n > no there holds 

s n (vo,@)<s n (@)<C/n. (15) 

Proof. Let vq G V W (S>) be any white vertex. If n > no(S>) is big enough, for each 
Z 2 -sublattice f2(vo) the set contains a [B(S)n\ -environment of vq in f2(vo), 



15 



where the constant B{S>) depends only on the construction parameters of Here 
denotes the largest integer smaller than pfE. Therefore we can choose s n (&) to be 
the maximum of S|_b(@)«j (^(v )) for all possible regular rhombic embeddings ^( vp) 



corresponding to Z -sublattice f2(vg). Now the claim follows from Corollary 3.3| for 
S'G-circle patterns. □ 

6 Uniqueness of quasicrystallic Z 7 -circle patterns 

In this section we consider quasicrystallic Z 7 -circle patterns as defined in [S] and then 
prove their rigidity. The proofs are based on the results of the previous sections and on 
similar arguments as for orthogonal Z 7 -circle patterns. 

6.1 Definition and useful properties 

Let tp € (0, 7r) be a fixed angle. Recall the definition of the labelling in We 
consider the following generalization of Definition |4.2| 



Definition 6.1 For < 7 < 2, the discrete map Z 1 : 1? + — > C is the solution of 

Q(fn,m: fji+l.m: fn+l,rn+li fn,m+l) ■ 

{fn,m ~ f n+1, m)(f n+1, m+1 ~ fn,m+l) ^2i(lj)— ir) (16) 

\fn+l,m fn+l,m+l*)(fn,m+l fn,m) 

and p| with the initial conditions 

Z 7 (0,0)=0, Z 7 (1,0) = 1, Z 7 (0,1) = e^-W. 

As in the orthogonal case one can again associate a circle pattern (for a quadrant 
of SG corresponding to 1? + and related to V and a^) to the map Z 1 , see [3] for more 
details. Furthermore, the following results generalize the orthogonal case. 

Theorem 6.2 ( 3J). (i) If R(z) denotes the radius function corresponding to the dis- 
crete conformal map Z 1 for some < 7 < 2, then 

(7 - l)(R(z) 2 - R(z - i)R(z + 1) - costPR{z)(R(z - i) - R(z + 1)) > (17) 

for allz£V\ {±N + iN\N G IN}. 



(ii) For < 7 < 2, the discrete conformal maps Z 1 given by Definition 4-2 are embed- 



ded. Consequently, the corresponding circle patterns are also embedded. 



Definition 6.1 can be generalized further. As explained in Section 13 of [SJ, discrete 
analogs of the power function z 1 can also be defined for quasicrystallic rhombic em- 
beddings $ instead of 1? + (or Z 2 ). In particular, let A — {±ai, ■ ■ ■ , ia^} C S 1 be the 
set of edge directions. Suppose that d > 1 and that any two non-opposite elements of 
A are linearly independent over R. For < 7 < 2 define the following values of the 
comparison function w on the coordinate semi-axis of lA: 

'l ifn = 0, 

w(ne k ) = { a 1 " 1 = e^" 1 ) 10 * * if n is odd, (18) 
TTH/-1 — — —t 2 - if n > 2 and n is even. 

The value of the logarithm loga^ is chosen as follows. Without loss of generality, we 
assume a circular order of the points of A on the positively oriented unit circle S 1 is 
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Figure 6: A quasicrystallic rhombic embedding with five-fold rotation symmetry. 





Figure 7: Example of a quasicrystallic Z 5 / 6 -circle pattern from the rhombic embedding 
of Figure [6j Construction (left) and corresponding circle pattern (right). 

a\, . . . , a^, — ffli, • ■ • , —a-d- Set a>k+d — ~ a k f° r k — 1, . . . , d and define a m for all m € Z by 
2d-periodicity. To each a m = e l9m € S 1 assign a certain value of the argument 9 m € R: 
choose 0i arbitrarily and then use the rule 

6 m +x - 9 m € (0, 7r) for all m e Z. 

Clearly we then have 9 m+ d = 9 m + ir. The points a m supplied with the arguments 9 m 
can be considered as belonging to the Riemann surface of the logarithmic function (i.e. 
a branched covering of the complex plane). Now, the branch of the logarithm is chosen 
such that 

log(ai) G [i9 m ,i6 m+d -i}, l = m,...,m + d-l. 

Using the Hirota Equation Q, this function w can be extended to the whole sector 
lA. Using suitable branches of the logarithm, w may also be extended to other sectors 
or to a branched covering of ifi . 

Figure [7] shows an example of such a quasicrystallic Z 7 -circle pattern. 

Note that for d — 2 the boundary conditions given in (|18|) lead to the circle patterns 



specified in Definition 6.1 Therefore, by Theorem 6.2 lul, the circle pattens corre- 
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Figure 8: A quasicrystallic rhombic embedding with ten-fold rotation symmetry. 





Figure 9: Example of a quasicrystallic Z 5 / 4 -circle pattern from the rhombic embedding 
of Figure [6j Construction (left) and corresponding circle pattern (right). 



sponding to the restriction of w to quad-surfaces C which are spanned by two 
coordinate semi-axis are embedded. We can apply finite and infinite flips to obtain 
other monotone quad-surfaces corresponding to rhombic embeddings. In particular, we 



5.4 



obtain restrictions to lA of the plane based quad-surfaces constructed in Example 
and Remark |5.6| Lemmas |5.5| and |5.7| imply that these lead again to embedded circle 
patterns. Thus we have proven 

Theorem 6.3 (Embeddedness of quasicrystallic Z 7 -circle patterns). Let il C be a 
simply connected monotone quad-surface. Then the circle pattern given by the function 



w with initial values ( 18 ) is embedded 



Example 6.4 (Construction of the examples in Figures [7) [9j and [IT]). The pictures 
in Figures [7J [9] and [TTJ have been obtained using a computer program implemented by 
Veronika Schreiber for her master thesis |24j . 

The rhombic embeddings in Figures [6] and [8] can be obtained by the construction 
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Figure 10: A quasicrystallic rhombic embedding with seven-fold rotation symmetry. 



method explained in Example |5.4| using the afBne planes E\ — {x 



:i + X\U\ + 

A2W2, Ai,A2 € H} and E2 — {x = ti + X1U1 + X2U2, Xi, X2 € H} respectively, where 
Ui, «2, t±, ti are defined as follows. {u±, U2, U3, U4, U5} is an orthonormal basis such that 



the matrix 




takes the form 



/cos(2tt/5) 
sin(27r/5) 



V 



- sin(27r/5) 
cos(27r/5) 






0\ 



cos(47r/5) -sin(47r/5) 

sin(47r/5) cos(47r/5) 

1/ 



with respect to this basis. The translation vectors are t\ = ( -0.2 -0.2 -0.2 -0.2 -0.2 ) 
and t% = (-0.5 -0.5 -0.5 -0.5 -0.5) respectively. Here v T is the transpose of the vector 
v. Translations in direction of the vector (11111) generally lead to rotationally 
symmetric rhombic embeddings. 

For the remaining construction it is important that the rhombic embeddings have 
rotational symmetry. Now consider the two sectors indicated by lines in Figures [6] and [H] 
or more precisely the corresponding 5-dimensional octant in 1? . Choose an exponent 
7 of the power function z 7 . In order to construct an image circle pattern which closes 
up as in Figures [7] (right) and [9] (right) we need to choose an integer p > 3 and take 
7 = 5/p. Given these ingredients, define the values of the comparison function w on the 
coordinate axes of the octant in Z 5 according to ( 18 1 and calculate the missing values 
for the octant using the Hirota equation Q. Taking the values on the quad-surface 
corresponding to the original rhombic embedding, we can construct a sector of the 
desired quasicrystallic circle pattern, see Figures [7] (left) and [9] (left). The closed circle 
patterns in Figures [7] (right) and [9] (right) are obtained using the rotational symmetry. 

The construction of the examples in Figures [10] and [Tl] is similar using a suitable 
affine plane in 17 . 
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Figure 11: Example of a quasicrystallic Z 7 / 5 -circle pattern from the rhombic embed- 
ding of Figure [6j Construction (left) and corresponding circle pattern (right). 

6.2 Uniqueness of quasicrystallic Z 7 -circle patterns 

In this section we prove uniqueness of quasicrystallic Z 7 -circle patterns using the same 
arguments as for the orthogonal case. 



We begin with a generalization of Proposition 4.6 Note that we need a (geometric) 
restriction which ensures that the kites corresponding to intersecting circles are convex. 
Unfortunately, flips (finite or infinite) may destroy convexity. 

Proposition 6.5. Let S> be a quasicrystallic rhombic embedding with associated graph 
G. Let a be the labelling corresponding to 3> . Let vq be an interior vertex of G with 
incident vertices V\ , . . . , v m . Consider two circle patterns for G and a with radius func- 
tions r and p. Denote Vj — r(vj) and pj — p(vj) for j — 0,1, ... , m and suppose that 
fj > rocosfx, and pj > pocosaij for j — l,...,m, where otj = a([vo,Vj]). Then 

3 = 1 HJ 3=1 PU 3 = 1 3 3 = 1 U 

where Cj = sin a j /((pj / po) + (po/pj) — 2 cos ay) for j = 1, ... ,m. 

Proof. The proof is based on a Taylor expansion of f a (x + logy) about y = 1. 

j- i , i ^ m i \ , hi \, t\ sin a(e x+los 5 - cos a) 2 

f a (x + logy = f a (x) + f' a (x)(y - 1 - ) '—^(y - if 

24^(cosh(x + log£) — cosa)^ 



by Lemma 2.3 (i) with £ = t + (1 — t)y for some suitable t € (0, 1) depending on x and 
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Equation Q for the two circle patterns implies 



n = zZ f°i ( log ^ = ^<7=i f<*t ( L 



Pa 



log « ) 



f r jPo 



- 1 



cos «j ) sin aj 



^ 2q(cosh(pj/p + log^j) - cosaj) 5 



where = + (1 



tj)^ > with suitable tj £ (0, 1) for j 



C _ i 

. . . , m. Furthermore 



cos a, = i 7 — + (1 - t,-) — 
Po r 



cos aj > 0. 



by our assumption. Thus 



rjPo 
r-oPj 



1 > 0. 



This implies the first claim since /^.(log = Cj by Lemma 2.3 (i) 



The second claim follows from the fact, that l/p and 1/r are also radius functions of 
circle patterns for G and a by Lemma 2.4 Also, the coefficients Cj are invariant under 
the transformation p i— > l/p. □ 

The following lemma specifies for which parameters the kites of the Z 7 -circle patterns 
are convex. For the cases excluded below, there exist non-convex kites, because already 
the kite built from the circle centered at the origin is non-convex. 

Lemma 6.6. If ip > tt/2 and 0<7<2 or if ip < tt/2 and (tt - 2ip)/(ir - ip) < 7 < 
tt /(tt — ip), all kites in the Z 1 -circle 'pattern given by Definition \6.1\ are convex. 

The proof is technical. A brief version is presented in the appendix [A] More details 
can be found in [T2"] . 

Theorem 6.7 (Rigidity of Z 7 -circle patterns from Definition [7|. If ip > tt/2 and < 
7 < 2 or if ip < tt/2 and (tt — 2ip) / (tt — ip) < 7 < tt / (tt — ip) then the Z 1 -circle pattern 
given by Definition 



6.1 



is the unique embedded SG -circle pattern for r Ir Jr and (up to 
global scaling) with the following properties. 

(i) The infinite sector {z = pe 1 ^ £ C : p > 0, f3 £ [0, 7(7r — ip)]} with angle j(tt — ip) 
is covered by the union of the corresponding kites of the circle pattern. 

(ii) The centers of the boundary circles lie on the boundary half lines. 

(Hi) All kites corresponding to intersecting circles are convex. 

The proof is a straight forward adaption of the proof of Theorem 



case) using Proposition 6.5 and Lemma 6.6 



4.4 



(orthogonal 



Theorem 6.8. Let ip £ (0, tt) and 7 £ (0,2) n , -E5]. Define Z^ -circle patterns 

on all four sectors Z± x Z± according to Definition 6.1 and glue these patterns to a 
circle pattern < &y on a cone with cone angle 2tt^. Then any embedded circle pattern with 
the same combinatorics and intersection angles which covers the same cone with one 
center of circle placed at the apex and which has only convex kites coincides with %? 7 ( up 
to scaling and rotation about the apex of the cone). 
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The proof is a simple generalization of the proof of Theorem 6.7 to circle patterns 
on a cone and combinatorics of Z 2 instead of 1?+. 

The following theorem is a direct consequence of the previous theorem and Lem- 
mas [5T5] and [5]7] on local changes of quasicrystallic circle patterns. 



Theorem 6.9 (Rigidity of quasicrystallic Z 7 -circle patterns I). Let S be a quasicrystal- 
lic rhombic embedding of a b- quad- graph covering the whole plane. Let A = {±ai, • ■ • ,±cid} C 
S 1 be the edge directions, where d > 1 and any two non-opposite elements of A are lin- 
early independent overR. Denote by ipmin the minimum of the undirected angles between 
any two elements of A. Let 7 € (0,2) with (tt — 2ip m i n )/(w — VWn) < 7 < tt / '(tt — tp m in) ■ 
Assume that the origin is a white vertex of 3>. Then a quasicrystallic Z 1 -circle pattern 
corresponding to S and embedded on a cone with cone angle 2ttj can be defined using 
the definition of the comparison function w on the 2d sectors of X d which contain the 



quad-surface f2@; see (18) and the remarks below. Assume further that the brick 11(0®) 
contains the whole lattice X d . 

Let be an embedded circle pattern with the same combinatorics and the same 
intersection angles which covers the same cone with one center of circle placed at the 
apex. Extend the comparison function w for If from fig to TL d . For each 1? -sublattice 
which contains two coordinate axes suppose that the corresponding circle pattern built 
according to this comparison function has only convex kites. Then c <a coincides with 
up to scaling and rotation about the apex of the cone. 

Note that the assumption on the convexity of the kites is only a restriction for a 



(small) neighborhood of the origin. This is due to Lemma 3.3 which implies that the 
ratio of the radii is almost one and thus the corresponding angles are almost the same 
in the isoradial case if the combinatorial distance to the origin is big enough. 

If all intersection angles of the labelling a are larger than tt/2, then the kites of 
any corresponding circle pattern are convex. Examples of such rhombic embeddings are 
suitable regular hexagonal patterns. Hexagonal circle patterns and in particular analogs 
of the holomorphic mappings z 1 have been studied by Bobenko and Hoffmann in [5] . 

Theorem 6.10 (Rigidity of quasicrystallic Z 7 -circle patterns II). Let & be a quasicrys- 
tallic rhombic embedding of a b-quad-graph. Assume that the corresponding labelling 
a : F{&) — > [tt/2, tt) has only values larger than Tt/2. Assume further that the origin is 
a white vertex. Let 7 € (0, 2). 

Define a quasicrystallic Z 1 -circle pattern for S> and a which is embedded on a 
cone with cone angle 2ttj using the definition of the comparison function on the 2d 



sectors of 1i d whose union contains the quad-surface see (18 1. In particular, the 
circle corresponding to the origin is centered at the apex of the cone. 

Let be an embedded circle pattern for Si and a which covers the same cone. Suppose 
that the circle corresponding to the origin is centered at the apex and that has only 
convex kites. Then coincides with ^ 7 up to scaling and rotation about the apex of the 



In this case that all intersection angles are larger than tt/2 the proof of Theorem 4.4 
can be directly adapted using Proposition |6.5| and the following lemma. 



Lemma 6.11. Let & be a quasicrystallic rhombic embedding of a b-quad-graph which 
covers the whole plane C Let G be the associated infinite graph built from white vertices. 
Then the simple random walk on G is recurrent. 

Proof. Recall that for the simple random walk all edges have conductance c(e) = 1 and 
also resistance r(e) = l/c(e) = 1. Our aim is to prove that the effective resistance R e g 
of the network G with these unit conductances is infinite. If two incident vertices of G 
are identified, that is the conductance of the connecting edge is increased to 00 and the 
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Figure 12: An illustration of the network obtained from G by shortening. 



resistance decreased to 0, then the effective resistance of the new network is certainly 
smaller. More generally, the proceedure of identifying a set of vertices of G is called 
shortening and reduces the effective resistance; see Q31 Section 2.2.2] or [551 Theorem 
(2.19)]. Therefore, it is sufficient to show that we have infinite effective resistance 
R' c g = oo for a network G' which is obtained from G by shortening. 

Without loss of generality we assume that all edges of @l have length one. Since *2i 
is quasicrystallic rhombic embedding the areas of the rhombi are uniformly bounded. 
Denote the lower bound by C\ > 0. Let vo € V(G) be any vertex and set Qk = 4fc 
for k E IN. Denote by V k C V(G) the set of vertices which are contained in the 
annulus A k = {z G C : g k -x < \z - v \ < g k } for k > 1. Then V(G) = Uf =1 V k . 
Identify the vertices of each V k to one new vertex v k - Then by construction, v k is 
only incident to Ufc_i for k > 2 and to v k+1 for k > 1; see Figure 12 Denote by \E k \ 



the number of edges which are incident to v k and v k +i for k > 1. Then the effective 
resistance of the shortened network is R' cS = YlkLi Furthermore \E k \ < area of 

{z € C : Qk — 2 < \z\ < Qk + 2}/Ci. Now a simple estimation shows that R' cS — oo. □ 
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A Appendix: Proof of Lemma 6.6 



First note that kites with intersection angle larger than 7r/2 are always convex. The 
remaining claims on convexity are consequences of Proposition 3 of [3] , namely 

(N+M)(R(z) 2 -R(z+l)R(z-i)~cos ipR(z)(R(z-i)-R(z+l)))(R(z+i)+R(z+l)) 

+ (M-N)(R(z) 2 -R(z+i)R(z+l)~cosipR(z)(R(z+i)-R(z + l)))(R(z + l)+R(z-i))=n (20) 



for ze¥ \ {N + iN\N e IN}, and the results of Theorem 6.2 



In particular, for tt — ip < tt/2 the kites with intersection angle ip < 7r /2 at black 



vertices are always convex. This follows from equation ( 20 ) and inequality ( 17 ) by simple 
calculations. This shows the claim for i/j > tt/2. 



If ip < tt/2 and (tt — 2ip) /(tt — xp) < 7 < it /(tt — ip), inequality ( 17 1 only implies that 
for all kites with white vertices z and z + i and intersection angle ip the angle at the 
point corresponding to z for < 7 < 1 and to z + i for 1 < 7 < 2 respectively is smaller 
than tt. This excludes some types of non-convex kites, but not all. 

For fixed n > 0, let T n be the piecewise linear curve formed by the segments 



[fn,m, fn,m+i] for m > 0. By Theorem 6.2 and its proof in [3], these curves are embedded 
without self-intersections and the vector v„(m) = f n ,m+i — fn.m rotates clockwise for 
< 7 < 1 and counterclockwise for 1 < 7 < 2 along T n . 
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Without loss of generality, we only consider the case 1 < 7 < 2 further. For < 7 < 1 
the proof is very similar. First, consider a kite on the symmetry axis, that is with white 
vertices corresponding to iK and i(K + 1). Then the assumption ^(n — ip) < it and 
the properties of the curves T n imply that the kites on the symmetry axis are convex. 
Next, we consider the intersection angles a and (3 of the line in direction of the vector 
v„(m) with the oriented half lines R + and e n ^~^^ 2 lR J + respectively. We additionally 
assume that n is odd. As the kites on the symmetry axis are convex and 7(7r — ip) < tt, 
we deduce a < ip and n/2 — ip < f3 < n/2. 

Finally, consider a kite with white vertex f n ^ m for odd n and m < n. We estimate the 
angles ct\ and a-i at this vertex of the kites containing the points f n . m -ii fn,m, fn+i.m 
and f n +i,m, fn,m, fn,m+i respectively. Note that these kites both have intersection 
angles ip. Using the above estimations on a and j3 and the properties of the curve 
r„ we obtain tt — 2ip < a\ and ai < tt. Therefore the angle at f n -i,m-i of the kite 
containing the points f n ,m-i, fn,m> /n+i,m is 2tt — 2ip ~ a.\ < tt. So this kite is convex. 
Furthermore, we deduce that the kite containing the points f n +i ym , f n ,mi fn,m+i is also 
convex. Consequently all kites containing the white vertex /„ !m are convex. 

Thus the Z 7 -circle pattern with ip < ir/2 and (tt — 2iP)/(tt — ip) < 7 < ir/(ir — ip) 
only contains convex kites. 
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